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The Effective Field Theory "without pions" at next-to-leading order is used to analyze universal 
bound state and scattering properties ol the 3- and 4-nucleon system. Results ol a variety ol phase 
shift equivalent nuclear potentials are presented for bound state properties of 3 H and 4 He, and for 
the singlet S-wave 3 He-neutron scattering length ao( 3 He-n). The calculations are performed with 
the Refined Resonating Group Method and include a full treatment of the Coulomb interaction and 
the leading-order 3-nucleon interaction. The results compare favorably with data and values from 
AV18(+UIX) model calculations. A new correlation between ao( 3 He-n) and the 3 H binding energy is 
found. Furthermore, we confirm at next-to-leading order the correlations, already found at leading- 
order, between the 3 H binding energy and the 3 H charge radius, and the Tjon line. With the 3 H 
binding energy as input, we get predictions of the Effective Field Theory "without pions" at next-to- 
leading order for the root mean square charge radius of 3 H of (1.6±0.2) fm, for the 4 He binding energy 
of (28 ±2.5) MeV, and for Re{a ( 3 He-n)} of (7.5 ± 0.6) fm. Including the Coulomb interaction, the 
splitting in binding energy between 3 H and 3 He is found to be (0.66 ± 0.03) MeV. The discrepancy 
to data of (0.10 =F 0.03) MeV is model independently attributed to higher order charge independence 
breaking interactions. We also demonstrate that different results for the same observable stem from 
higher order effects, and carefully assess that numerical uncertainties are negligible. Our results 
demonstrate the convergence and usefulness of the pion-less theory at next-to-leading order in the 
4 He channel. We conclude that no 4-nucleon interaction is needed to renormalize the theory at 
next-to-leading order in the 4-nucleon sector. 
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I. INTRODUCTION 

Nuclear processes, both with and without external probes, at energies well below thepion production threshold 
are naturally described by the Effective Field Theory "without pions" EFTy- (see e.g. [IJ-Q for reviews). It shares 
with QCD the symmetries and effective low-energy degrees of freedom and provides a simple, systematic, and model 
independent approach appropriate for systems at very low energies. Like all Effective Field Theories, it allows for a 
systematic expansion of the scattering amplitude in a small, dimensionless parameter to attain the desired level of 
accuracy in observables and thereforeprovides reliable a priori error estimates. Its simple structure makes it also 
easier to handle than pionful EFT jg-|l|. 

EFTy. has successfully been applied to many reactions of two and three nucleons with and without electro-weak 
currents. While calculations at fourth order have been employed to yield very accurate results in those systems (see 
e -9- @)i higher order calculations involving more than three nucleons are still in their infancy. This work supplements 
the few, already existing leading-order (LO) calculations in the bound four-nucleon system [ToL |Tl[ by a next-to- 
leading order (NLO) analysis including the Coulomb interaction and the leading-order 3-nucleon interaction. As 
first observable in the four-nucleon scattering system, it adds the singlet S-wave 3-Helium-neutron scattering length 
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ao( 3 He-n). By that, it addresses the question whether EFT^f still converges for the a-particle. 

Our investigation has three central goals: (i) to support evidence that the four-nucleon bound- and scattering system 
can be described by EFT^ ; (ii) to provide strong evidence that no four-nucleon- interaction is necessary to renormalize 
EFT^ at NLO; and most importantly (iii) to clearly demonstrate the feasibility of the variational approach of the 
"RRGM" (RRGM) for EFT^ in addressing these goals also for future studies involving heavier systems. 

EFT^f systematically expands observables in powers of the typical low-momentum scale pt yp of the process consid- 
ered, measured in units of the scale Aeft^ at which the pion can be resolved as a dynamical degree of freedom, set 
by the pion mass m n . Equivalently, the expansion parameter is given by the ratio of the range w 1.4 fm of the 
one-pion-exchange over typical resolutions or sizes of the process under consideration. The 4-Helium binding energy 
is 28.5 MeV and its typical size of about 1.6 fm is comparable with the approximate range of the one-pion-exchange. 
A s yst ematic expansion in the four-nucleon system could therefore converge very slowly at best. However, Platter et 
al. [10| showed that EFTjf at LO can explain with good accuracy the correlation between the 4-Helium and triton 
binding energies, the famous Tjon line [12|]. One focus of this article is to test whether or not EFT^ at NLO shows a 
convergence order-by-order, and converges to experimental values. 

For this purpose, we investigate correlations between three- and four-nucleon observables. They correspond to 
universality classes of nuclear interaction models, all resulting in on-shell (phase shift) equivalent 2-nucleon potentials 
which share the same two-body scattering lengths a s> t in the singlet and triplet S'-wave. The elements of a class, 
parameterized by a s _ t , therefore fully describe the nucleon-nucleon (NN) system in the zero energy limit. In con- 
tradistinction, nonuniversal parameters like the effective range only add perturbative corrections and correspond to 
higher order interactions in the EFT pow er-counting. Their contribution to an observable vanishes in the zero energy 
limit. The empirically found Phillips [13[ line, the Efimov spectrum of three-body bound states and the Tjon line [12[ 
are examples for universal properties of NN potentials of the same universality class defined by a Sj t (for reviews see 
e.g. QJ]). The Phillips line is a correlation between the triton binding energy B(t) and the neutron-deuteron 

doublet scattering length and can be explained as a result of the fact that on-shell equivalent 2-nucleon potentials 
have in general different off-shell behavior. Such a dependence of low-energy observables on details of short-distance 
2-nucleon physics is however unphysical. This renormalization defect can be corrected by including a three-body 
interaction which eliminates the dependence on the NN potential. The datum on the Phillips line is then reproduced 
at LO by fitting the strength of the three-body interaction [l6[. In the zero energy limit, the 3-nucleon system is 
hence described by three universal parameters: a Sl t and one three-body low-energy datum. Different values for this 
low-energy datum result in a different strength of the 3-nucleon interaction. As the latter is however also dependent 
on the NN potential used, one can use an alternative approach to explore the dependence of few-nucleon observables 
on the additional three-body input: Set the three-body interaction to zero and map out the dependence of observ- 
ables on the regularization procedure used for the 2-nucleon potential. The results are correlations between few-body 
observables. In this work, we employ all three procedures to map out correlations: different values of the regulating 
cutoff, a variation of unresolved short-distance observables by allowing for non-zero, higher order P-wave and SD 
interactions, and a variation of the coupling strength of the three-body counterterm with fixed 2-nucleon interaction. 

One might be tempted to choose form and value of the cutoff for the 2-nucleon potential such that the triton binding 
energy is reproduced "exactly", and then perform all calculations with such a highly fine-tuned potential, see e.g. [l7| . 
However, such an approach would deprive one of a simple procedure to find an estimate of the theoretical uncertainty 
of the calculation. Including electromagnetic interactions will also be more cumbersome, as these form-factors have 
to be gauged and additional, transverse currents need to be constructed. On top of that, it is not clear that one can 
find such a 2-nucleon potential to eliminate any 3-nucleon interaction, or that its choice would be unique [64|. The 
numerical advantage of not having to deal with 3-nucleon interactions in some systems can in general be balanced by 
the problem to find the "right" 2-nucleon potential. In any case, such additional effort does not increase the reliability 
of the calculation. We therefore keep the cutoff arbitrary to map out correlations and error estimates. 

Adding one more nucleon, the Tjon line of a one-parameter correlation between the binding energies of the triton 
and 4-Helium supports the assertion that no additional universal parameter is necessary to classify the four-nucleon 
system. Its emergence in EFTjf at LO was demonstrated by Platter et al. [10]. In order to extend this finding to 
NLO, this work investigates not only the dependence of the triton charge radius, the splitting in the trinucleon binding 
energies, and the 4-Helium binding energy on the triton binding energy, but also finds a similar connection between 
the real part of ao( 3 He-n) and the triton binding energy. Therefore, it supports the expectation, based on naive 
dimensional analysis, that no four-body interaction is necessary to renormalize the theory at NLO in the four-body 
sector. Up to this level of accuracy, observables in the four-nucleon system are universal consequences of the two- and 
three-body dynamics of its constituents. For all observables we consider, namely the charge radius of the triton, the 
trinucleon binding energy splitting, the 4-Helium binding energy B(a) and the scattering length ao( 3 He-n), we find 
predictions consistent with experiment within error margins when the 3-nucleon interaction strength is fitted to B(t). 

For this work, the Refined Resonating Group Method [l8| is employed. This variational method numerically 
solves the Schrodinger equation in coordinate representation. Within a model space spanned by Gaussian type 
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wave functions, minimization of a respective functional yields scattering and bound state observables for a potential 
parameterized in terms of Gaussians. With the RRGM, speedy scattering calculations are feasible in the A — A system 
and even beyond, because of the relatively simple structure of the NLO EFTjf potentials. To study the feasibility of 
using EFT^ in the RRGM is the other focus of this article. In the course of this investigation, it is also imperative to 
carefully assess that the spread of results does not come from numerical inaccuracies. 

The article is organized as follows: First, a brief overview of EFTjf is given as a low-energy theory which allows 
for a systematic improvement of the interaction amongst nucleons up to the desired accuracy. An introduction to the 
calculational tool of the RRGM in sect. Mil is followed by the derivation, structure, and fit of the NN potential in 
sect. II VI The results are presented in sect. El followed by a concluding section. An appendix addresses the numerical 
accuracy and computational costs of the calculations. 



II. PION-LESS THEORY 



In this section, the theoretical framework of an Effective Field Theory for nucleons without pions is recapitulated. 
For details, the reader is referred to the exemplary reviews [!, Hp. 

The nuclear potential is derived from a Lagrangean (see e.g. |20||) with contact interactions which are momentum 
independent at LO and momentum dependent at NLO: 



-f,NN 
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The only effective low-energy degree of freedom is the nonrelativistic fermion iso-doublet N = °f Weyl spinors 

p and n for the proton and neutron, respectively. In eqs. (|1I2[) . Einstein's summation convention is understood, i.e., 
a sum from 1 to 3 over repeated Arabic indices. The Lagrangean is an isoscalar at this order, with neutrons and 
protons having the same mass M. The coupling constants c| N ) LO are referred to as low-energy constants (LECs) or 
counterterms, and Oi arc the Pauli spin matrices. The power-counting in EFT^ at NLO (see e.g. results in four 
independent LECs, with zero P-wave and SD-transition amplitudes providing the five constraining equations (see the 
discussion of eq. (jTTJ) for how this constraint is implemented) . 

Within the EFT framework, it was found [ml2l| - l25j that exactly one 3-nucleon contact interaction (3NI) is necessary 
to renormalize the A = 3 system in the doublet-S channel at LO and NLO. In our NLO calculation, this three-body 
counterterm has the form 

C f , 3 N = C^ 1 (N^N)(N^nN)(N^T l N) . (2) 

Other forms are identical after Fiertz-transformations [22]. Once the 2-nucleon parameters in eq. ((T|) are fixed, the 
universal correlation lines between three- and four-nucleon observables, discussed in sect.fVl are parameterized by the 
3-nucleon coupling CJnj. 

The amplitudes derived from the Lagrangean in eq. ([!]) are given as an expansion in the dimensionless parameter 
Q ~ Ptyp/Ab, where pt yp is a typical low-momentum scale of the system, and a rough estimate for the breakdown scale 
Ab of the theory is the pion mass. In the neutron-proton system at center of mass energies of less than the deuteron 
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binding energy -B(d), pt yp is set by the binding momentum, \J M B(d), which leads to an expansion parameter of 
Q < |. A calculation at N"LO is then expected to be accurate up to perturbative corrections of order Q n+1 . For 
energies below B(d), the parameter Q is found to be approximately constant (e.g. [26j]) but increases for higher 
energies. This increase leads eventually to a breakdown of the perturbative expansion. In practice, a cornucopia of 
EFT^ calculations 0, H, HH H3, HH to higher orders has shown that the expansion converges somewhat beyond the 
pion mass and that i < Q < |. 

The NN interaction is understood as an effective potential following Weinberg's original definition [2{|. In the 
pion-less theory, the potential reduces to tree diagrams of two in- or out-going nucleons, with vertices from eq. (fT|). 
While not necessary from the EFT standpoint, it is convenient for this work to insert the full NLO potential in this 
form into the Schrodinger equation instead of treating the non-leading terms in perturbation theory. This course of 
action has been pursued regularly and includes some contributions which are formally of higher order in the EFTjf 
power-counting but does not increase the accuracy of the result. The re-summation does however lead to spurious 
bound states in the two-particle sector which can impact three-particle observables, see e.g. [28| . 

This method also allows us to briefly comment on how our results would change in an alternative power-counting in 
EFTjf proposed by Beane and Savage [3(|, where both scattering lengths and effective ranges count as l/pt yp ~ Q" 1 ■ 
In that case, effective-range corrections must be re-summed, and some combinations of the NN interactions 
in eq. ([T]) are promoted to LO. One would now be compelled to iterate the potential in the Schrodinger equation. 
As one performs technically the same steps as above, one arrives at the same amplitude. Our results can thus also 
be interpreted as LO calculations in the alternative formulation. The difference between this alternative and the 
approach taken here is therefore only in the question whether re-summing effective range contributions is optional or 
mandatory. If the effective ranges can but do not have to be included as NLO corrections, one expects that observables 
change only by parametrically small amounts when one calculates first with zero effective ranges (i.e. keeping C^® 
only), and then adds the supposed NLO terms Cf^. Our results for the triton charge radius in sect. IV Al and for 
the Tjon line in sect. IV CI confirm this assumption in the three- and four-nucleon system. From that perspective, we 
therefore see no reason to make the effective-range resummation mandatory. 

As pointed out in [3l| . for a short ranged potential like this, the Wigner bound can potentially constrain the value 
of the effective ranges: When the scattering length is positive, an energy-independent potential exists only below a 
cutoff-dependent upper bound for the effective range. For high enough cutoffs, the physically observed effective range 
in the 3 Si channel will exceed this Wigner bound. However, we will demonstrate in sect. IIVI that for our choice of 
cutoffs, the problem does not arise. 

The derivation of the explicit form of the potential in coordinate space from eq. ([T|) is postponed to sect. IIVI in favor 
of first introducing the numerical method used, as it motivates also the choice of operator structure for the potential. 



III. RESONATING GROUP METHOD 

We employ the variational method of the Refined Resonating Group Method (RRGM) [l8[ to solve the Schrodinger 
equation in coordinate space. In this section, the method is introduced. Factors determining the numerical stability 
of EFT,f calculations with the RRGM are addressed in app. [A] 

The RRGM uses a Gaussian basis to span the variational space and a Gaussian expansion of the radial dependences 
of the nuclear potential which allow for an analytic calculation of the Hamilton matrix. As in every variational 
approach, the basis is incomplete. Care has to be taken to avoid linearly dependent basis vectors. The computer 
time for the calculation depends on the dimension of the variational space, the number of Gaussians needed for an 
accurate fit of the radial functions of the potential, and its operator structure. A Gaussian regulator for the contact 
interactions leads directly to a Gaussian radial dependence of the potential and makes the RRGM an efficient tool to 
analyze the few-nucleon sector. This section describes those aspects of the method which are necessary to understand 
the choices made for the form of the regulator function, for the range of its cutoff values, and for the set of operators 
which constitute the potential. 

To determine the bound state wave function, the Ritz functional is minimized in a model space spanned by vectors 
of the form 




N-l 

[n 

L fc=l 



i J 



e- likP *y hi {p k ) 



(3) 



The system consists of N particles. To each of the N — 1 Jacobi coordinates pk, one assigns a set of different 
width parameters jdk and of different angular momenta Iki represented by solid spherical harmonics [32]. The 
antisymmctrizer is denoted by A, and the square brackets indicate angular momentum couplings, with the orbital 
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part Li being combined with the spin Sj to the total angular momentum J. The label d distinguishes different sets 
of width parameters 7, while i and j label sets of orbital- and spin angular momentum coupling schemes. The spin 
function S is constructed as a sum of products of the single particle spin functions s. The isospin function T is built 
analogously from single particle isospin functions and distinguishes between neutrons and protons. The superposition 
coefficients Cdij are determined by minimizing the Ritz functional. The magnetic quantum numbers of the spherical 
tensors are not explicit in the above equation where they are not needed to label a specific basis vector. Since one 
deals only with reduced matrix elements, even the magnetic quantum number corresponding to the total angular 
momentum J is of no significance. 

For the scattering state, the Kohn-Hulthen variational principle is used with the following ansatz for the wave 
function: 



= A< 



E 



7*1 j 



SxjF Lj (Rj) + a Xj G Lj (Rj) + £ bx jm R} ,+1 e- 



(4) 



J 1,a 

The two fragments are represented by bound state wave functions r/> • 1,2 , determined by the aforementioned Ritz 

minimization and built from vectors as in eq. The scattered fragments are separated by R, and the orbital 

angular momentum between the two fragments is carried by a spherical harmonic Yj,. A specifies the boundary 

condition which allows regular Coulomb waves Fl [33| only for the channel j = A. The reactance coefficients a\j 

and the b\j m are determined by minimizing the Kohn-Hulthen functional, taking into account rif. channels. The 

additional set of variational parameters b\j m is necessary to approximate the wave function in the interaction region. 

The irregular Coulomb functions are regularized with a polynomial weighted by an exponential and renamed Gl- 

To become more familiar with the terminology and as a precursor to the calculation presented in sect. El consider 

a neutron scattered off a 3-Helium nucleus. In that case, rik = 654 channels were included, the fragment wave 

1+ 1 + 

function of the neutron is ip?^ = 1, and that of the 3-Helium ip^^y ^he letter consisted of 224 terms, namely 78 

(Li = 0, L2 = 0) components, called SS— configuration, 82 SD-, 45 DD-, and 21 PP— configurations. An example 
of an SD-configuration might be given by the following parameters: The angular momenta l\ = 0, 12 = 2 on the two 
Jacobi coordinates couple to a total orbital angular momentum L = 2. This mandates the individual spins of the 
nucleons to be aligned, i.e., total spin S = §, so that finally a part of the triton system with JJ = i + is formed. 
The two total angular momenta of the fragments, Jf 2 = j + , can be coupled to channel spins, S c 6 {0, 1}. However, 
only S c = has to be included for the — + -channel, which in turn dictates that the orbital angular momentum 
between the fragments is zero, Lj — 0. The width parameter u>j m for the inter-fragment wavefunction plus two widths 
7dfe for the 3-Helium bound state component complete the specification of the basis vector. In the following, the 
model space is defined to be a vector space spanned by basis vectors V'bs/SS as given in eq. or eq. (UJ). 

Once the potential is written in terms of spherical tensor operators and its radial dependences are expressed in 
Gaussian functions, all the coordinate space matrix elements which are needed to minimize the respective functional 
can be cast into the form 



/ = / d 3 Pl . . 



■d 3 Pn k - 



1 exp 



rik — 1 



ni 

- J2 c» v 'p v -p v ,+Y.?-pv 



(5) 



The matrix C uv transforms the Jacobi coordinates in the in-coming channel and the relative coordinates of the radial 
dependences in the potential operator to the coordinates in the out-going channel. Furthermore, it takes into account 
that the antisymmetrizer permutes single particle coordinates. The vector is related to the generating function of 
spherical harmonics. Its dimension is n/, namely the number of spherical harmonics in the matrix element /. The 
integrals I can hence be evaluated analytically for all operators of the EFTjf potentials and the Coulomb interaction. 

For an accurate description of the scattering state, so-called distortion channels have to be added to the physical 
channels in eq. ^ . They increase the variational space to allow for a more accurate description of the wave function 
in the interaction region. There, a separation into two bound fragments as mimicked by the physical channels does not 
resemble, e.g. three- or four-body breakup states. Those channels do not have an asymptotic tail, i.e., only the square 
integrable terms with coefficients b\j m in eq. ([4]) constitute the relative part of their wave function. Furthermore, 
their fragment functions are not restricted to describe bound states, but only to have the correct quantum numbers. 
That means it suffices to specify a single set {Iki, "fdk, Sj}, instead of superimposing multiple sets to a bound state. 
In the calculations, all but one of the components of each of the physical channels were used for this purpose, e.g., a 
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physical 3-Helium-neutron channel, with 224 basis vectors for the 3-Helium fragment, yields 223 configurations. The 
number of distortion channels, and by that the dimension of the model space, is ultimately determined by the number 
of different w jm used in eq. (j4|). With the number of physical channels fixed, more distortion channels were added to 
obtain a converged result for the observables. 

In view of the error analyses performed in sect.|Vl it is imperative to demonstrate that purely numerical inaccuracies 
of the variational method are no significant source of error. In app.EJ we demonstrate that differences in calculations 
of the same observable do indeed not originate from numerical inaccuracies. 



IV. TWO-NUCLEON POTENTIAL AND PARAMETER DETERMINATION 

In this section, the potential in coordinate representation is derived from the Lagrangean in eq. ([I]), following the 
lines of [13] but choosing a different operator structure which simplifies the implementation into the RRGM. In the 
second part, the fitting procedure of the low-energy constants (LEC) and the experimental input is explained. 

In momentum representation, the two-body potential for an A-nucleon system following from the nine four-nucleon 
contact interactions of eq. ([TJ is 



rANLO) ' 

EFTjf ,NN 

i<j 



E (cl° + c\° Si ■ ffj + cf L V + c 2 NLO fc 2 + Si ■ Si (cf L V + c 4 NLO £ 2 ) 

+ iC ,NLO (fi+M .q X k + C 6 NLO (f • SiS- Sj + C 7 NLO fc • Sik ■ S 3 } , (6) 



with q = p — jf , k = p ^ defined in terms of the in(out)going center of mass momenta p(j?) of one nucleon. 
Following eq. ([2]), the three-body potential is given as 

A 

^EFTy,3iV = E QjNI ( T « ' T j + T k ' Ti + Tj ■ ffe) . (7) 
i<j<k 

Regularized with f\(q) = cxp (— (^/A 2 ) and Fourier transformed, eq. ([5]) is cast into the form 



Io(A,r) (A 5 + A 6 ai ■3 j )r z + 

Jo (A, r) A 7 L -S + I (A, r) A 8 (ffi ■ ra 3 ■ r - ^r 2 ^ • B. 



o (A, r) , [ [cT ® s ] 2 ® [of ® erj] 2 ] °° i , (8) 



and eq. ([7]) into 



^EFT^ ,3JV = E 7 " ( A ' ^1 ) 7 " ( A ' ^ ) C 3NI Ti ■ T, , (9) 

i<j<k 

cyclic 

with the interparticle vector r — ri — fj, the orbital angular momentum operator L = —irx V, and total spin operator 
S = \ {S-i + S.j ■) , as well as the regulator function I Q (A,r) = exp (— A 2 r 2 /4). It is not necessary to symmetrize the 
potential with respect to particle exchange because the basis states are antisymmetric. Instead of fitting the original 
LECs Ci, it was favorable to adjust directly linear combinations 

9 

A i = Y,*ijV l{i) Cj ■ (10) 
j'=i 

This re-definition avoids a fine-tuning of large versus small terms in the sum which arises because the anticommutators 
absorb powers of A into the derivative acting on the regulator function I Q (A, r). 

The operator sets of AVf8 and CD-Bonn are subsets of the one of eq. ([5]), while the relatively complicated radial 
dependences of the former contrasts with the polynomial- weighted Gaussians of eq. ([5]) . Although most of the operator 
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structure is hidden in the anti-commutators, Hermitecity is manifest in this form, and the Gaussian radial dependences 
allow for efficient RRGM calculations. The natural size of the coefficients A t is not trivial to estimate because the 
derivatives in the two anticommutators act not only on the wave function but also on the Gaussians in the potential. 
Therefore, the parameters ^3,4,9 are of lower power in A and hence differ considerably in size relative to the others. 



can be employed to reduce the number of parameters to five. Here, we do not impose those constraints exactly, but 
instead extend the x 2 -measure to fit the P-wave phase shifts to a fraction of the Nijmegen values, 



This is fully consistent with the EFT philosophy that higher order interactions can only induce higher order corrections 
in observables. Different P-wavc constraints, compatible with eq. (fT2"|) . allow for such a controlled modification of 
short-distance structure. An exact implementation of eq. (|f f P provides no significant gain in computer time relative to 
eq. (fl"2"|) . We therefore exploit the additional handle on higher order effects provided by eq. (fT2|) to gauge the accuracy 
of a NLO calculation. A consequence of this approach are non- vanishing rank one and two interactions corresponding 
to LECs A7 and As. 9, respectively. These tensor structures are also found when one implements the SD-interactions, 
which enter beyond NLO. We choose to include these as representations of higher-order effects and to constrain their 
parameters by the SD-mixing angle. 

Different sets of weight factors for the phase- and -B(<i)-deviations were used for the various potentials. Therefore, 
this x 2 is n0 objective criterion for the quality of the potential and we abstain from quoting it. 

The basis of this potential does not mix the singlet with the triplet NN channels. Therefore, it was convenient to 
fit the projections of the potential in the spin singlet- and triplet channels separately. In the former, three parameters 
were adjusted to reproduce 5Nij( 1 S'o) and <5fit( 1 -Pi)- In the triplet channel, the fit is initially to B(d) only, followed by a 
fine-tuning of the six LECs to refine Sat( 3 Si), <5fit( 3 -fo,i,2) and the SD-mixing angle eifit( 3 Si— 3 Di). Here, the following 
hierarchy of weights u; (observable) in the \ 2 ~ function was used: w(B(d)) > w ((5( 3,1 5'i i o, £i)) > w (S( 2S+1 Pj)) , i.e. 
the most weight was put onto B(d), and we typically used a ratio of w(B(d))/w(5) of the order of 10. Dependence 
of the resulting values for the nine LECs on the model space was minimized by using an almost complete set of 40 
Gaussian basis states as described in sect. Mil 

The ^-minimization which determined the LECs was carried out with a modified version of the genetic search 
algorithm already used for wave function optimizations in [34j ] . This algorithm has the advantage of being independent 
from an educated guess for an initial set of parameters. With the terminology specified in [341 ], the search parameters 
for the algorithm were chosen as follows. The initial population was set to consist of more than 10 4 individuals, with 
each individuum corresponding to a set of LECs. This relatively large number should ensure that a good fraction of the 
entire parameter space is probed. The search intervals used for the Ai's were chosen differently: ^3^.9 s [— 10 3 ; +10 3 ], 
while ^1,2,5— 8 S [ — 10 4 ; +10 4 ]. A factor of up to A 4 is not included in the ^3,4,9 but contributes when the derivatives 
in the anticommutator act on the regulator functions (see above). Therefore, we expect the ^3,4,9 to be smaller than 
the other A^s. 

The contribution a specific operator made to B(d) differed amongst the l Vf. While some EFTjf potentials distribute 
B{d) similar to the potential models, AV18 [35| and CD-Bonn 36], i.e. the largest fractions coming from the central, 
and tensor terms, others summed B{d) quite differently. However, there is no physical reason why one operator of 
l V J f should contribute more than the other. Not only are LO and NLO operators combined in (r) (see eq. ©), 

but they also have portions of the two central operators hidden in the central anticommutator, as well as a tensor 
component in the rank two anticommutator. 

A heuristic explanation for the emergence of a variety of different LEC sets each representing a valid EFT^ inter- 
action shall now be given. The internal structure of the nucleons is encoded in the LECs, and different regulators 
lead, after proper renormalization, to the same low-energy NN observables. Therefore, no unique set of coupling 
constants is expected. Here, two methods are employed to model different short-distance interactions with the same 
long distance behavior. First, the regulator was altered by changing the magnitude of the cutoff. Second, different 
sets of LECs were found by using different input for the fit. 

The difference in input data must be compatible with the order at which the EFT calculation is carried out in the 
following sense. A consistent EFT calculation at order n predicts low-energy observables accurately up to uncertainties 
of order Q n+1 in the dimensionless expansion parameter. This property extents to the data used to fit the LEC, i.e. 
there is no need to refine the fit of the parameters in table U for a better reproduction of, e.g. the deuteron binding 
energy B(d) of the potentials. By weighting data differently in the x 2 -function, different sets of LECs were found, 
see table |U The values for B(d) and the neutron-proton phase shifts for those sets were all in the NLO uncertainty 




(11) 




(12) 
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TABLE I: Numerical values of the cutoff A and the LECs of the NLO EFTy potentials. The LECs were fitted to B exp (d) [37 
and NN phase shifts. 





A 


Ai 


A 2 


A 3 


A 4 


A 5 


A 6 


A 7 


A 8 


A 9 




[MeV] 


[MeV] 


[MeV] [MeV-fm 2 ] 


[MeV-fm 2 ] 


[MeV-fm -2 ] 


[MeV-fm- 2 ] 


[MeV] 


[MeV-fm" 2 ] 


[MeV-fm 2 ] 




414 


-143 


464 


-57.0 


68.8 


-52.4 


-346 


-3.97 


-62.1 


-0.111 


7* 


432 


-612 


944 


-157 


110 


351 


-723 


66.0 


-168 


-0.137 




544 


-1224 


1036 


-432 


336 


1704 


-1851 


-78.9 


-378 


-0.130 




544 


-1625 


-89.6 


-298 


49.8 


1870 


34.3 


-990 


-734 


2.13 ■ 10" 3 


7* 


544 


-405 


317 


-125 


108 


287 


-493 


-271 


-360 


-9.88 ■ 10~ 4 


7* 


648 


84.9 


-324 


-888 


404 


4845 


-2064 


-342 


-2319 


-0.716 




648 


-1316 


1039 


-530 


431 


3775 


-3666 


-75.4 


-1052 


-0.178 




672 


301 


-143 


-708 


480 


3527 


-3450 


-192 


-1613 


-0.257 


9 v, 


672 


-158 


-201 


-362 


275 


2001 


-1659 


79.7 


-1270 


-0.0850 



TABLE II: Two-nucleon observables calculated with NLO EFT^ potentials as defined in eq. (|8]) and table U The deuteron 
binding energies B(d) are labeled by the RRGM model space in which they were calculated. The data for the scattering lengths 
a Stt and the effective ranges r St t is taken from [3Sl ]. 



> > > > > > 

a) o> a) a) e a> 

s s s s s s 



exp - 2.225 2.225 2.225 -0.120 0.0275 0.0138 0.00892 

Wf 414 2.239 2.238 2.226 -0.121 0.0275 0.0136 0.00914 

2 V f 432 2.224 2.212 2.203 -0.119 0.0279 0.0134 0.0104 

544 2.221 2.205 2.153 -0.123 0.0268 0.0144 0.00824 

A V i 544 2.233 2.201 2.185 -0.119 0.0270 0.0133 0.00859 

5 Vf 544 2.224 2.222 2.186 -0.119 0.0258 0.0133 0.00666 

6 V ?t 648 2.212 2.019 2.034 -0.120 0.0271 0.0132 0.00835 

7 V f 648 2.214 2.137 2.012 -0.115 0.0260 0.0124 0.00687 

*Vf 672 2.221 1.985 2.124 -0.118 0.0264 0.0125 0.0077 

9 V # 672 2.250 2.182 2.114 -0.118 0.0255 0.0125 0.00648 



range around data, see table HU Another approach which yields potentials with the same long-distance but different 
short-distance structure is changing the input for the P-waves. The P-wave phase shifts are observables of higher 
order, and therefore the values predicted by the NLO potentials are only restricted to a range compatible with zero. 
We enforce different P-wave interactions by including and varying an appropriate term in the x 2 -function, see eq. (|12[) . 
We employed a combination of both methods, different weights and different P-wave input, to find several NN phase 
shift equivalent potentials for a single cutoff value. 

Another point had to be kept in mind for the fit. In principle, the EFT philosophy allows for deeply bound states 
because reactions at energies for which the pion-less theory is applicable will not probe those low-lying, nonphysical 
states. As mentioned in sect, [ill they can occur when re-summing effective range contributions. Nevertheless, we 
considered only potentials which were found not to posses any of these deeper lying bound states. The reason is that 
the variational basis does not span the entire Hilbert space and therefore might be insufficient to expand a low-lying 
and hence very localized ground state of a potential in the course of the fit. In that case, it would be impossible 
to tell if this state is low enough in energy to be considered marginal. To ensure that no such ghost states were 
present during the fit, vectors with narrower width parameters were added to the basis for the deuteron and the 
triton to model more localized wave functions and sort out the troublesome LEC parameter sets. Unphysical states 
might also be formed as clusters in a three- or four-nucleon calculation because of narrower basis states, in which 
case the corresponding LEC set was sorted out. We employed as figures of merit for a potential a converged value of 
the x 2- function after successive runs of the genetic algorithm, and the numerical stability of the deuteron and triton 
binding energies with respect to changes in the model space, see discussion below. The numerical values of the 
LECs over a range of cutoffs are presented in table [H The results for the deuteron binding energy in three different 
model spaces per potential are listed in table [TT] to demonstrate the approximate closure of the set u>i2o which was 
used to fit the LECs. A model space dependence of the potential would lead to different deuteron binding energies 
B(d) in W12 and 10120- Those sets were 40-dimensional, while the set wqs is only 9-dimensional, with 6(3) optimized 
widths for the S(D) wave function component. It was used in the A — A calculations to build the deuteron fragments. 
The change in B{d) associated with the substitution of u>i 2 o by w\i increases with the cutoff. This is due to the 
absence of narrower width parameters in w\2 which, as seen above, has a greater impact for narrower potentials, or 
larger cutoffs, respectively. Hence, this change is attributed to a shortcoming of the w\2 model space rather than to 
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one of the potential. For the potential s Vf, B(d) decreases considerably from W120 to 1012 but increases again and 
approaches the u>i2o value in a much smaller but optimized space wq^ . In general, the quality of the model space used 
for the fitting procedure is comparable to the one used for the four-nucleon calculations presented in the next section. 

The decrease in B(d) when going from W120 to wqs has no substantial impact on the results for two reasons. First, 
u>63 was only used to expand the deuteron cluster in the 4- Helium system and not the one in the triton/3- Helium, 
which was built and optimized independently of the deuteron parameters. Second, in the A = 4 system only scattering 
reactions at energies below the deuteron-deuteron threshold were considered. Therefore, the deuteron configurations 
contribute only through the distortion channels. This leads to the non-converged results for the imaginary part of the 
4-nucleon scattering length, see sect. IV Dl 

The width parameters which constitute set Wes were tailor-made for each potential using the genetic algorithm. 
In this process, the found for smaller cutoff values consisted on average of wider Gaussian width parameters 
compared to the ones for relatively large values of A. This behavior was already explained in sect. II VI in the discussion 
of fig. El While B(d) was used as experimental input to determine the LECs, the singlet and triplet scattering lengths 
a s .t were used only implicitly through their dependence on the phase shifts in the respective channels. They were 
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calculated from the appropriate S-matrix elements at a center of mass energy, E cm = 1 keV (see discussion below, 
eq. (IB)). " 

As an example, the phase shifts for x Vf are shown in fig.[T]along with the phases from the Nijmegen analysis |30L t40l| . 
Only phase shifts for a center of mass energy below 0.35 MeV were used to determine the LECs, so that phases for 
higher energies can be used to gauge the quality of the fit. The other potentials 2 ~ 9 Vjt reproduce the Nijmegen phases 
with the same accuracy, i.e., the phases up to 10 MeV differ by less than 1%. While <5Nij( 1 S'o) and 6mj( 3 Si) are 
reproduced up to 1% for E cm < 10 MeV, the SD-mixing angle e\ reaches, in the same energy interval, only a 10% 
accuracy and also showed a stronger variation between the various potentials. Both uncertainties are within the error 
margins predicted by EFT^, i.e., approximately 10% for the S- waves at typical momenta ptyp ~ 45 MeV and 30% for 
the higher-order observable t\. For all potentials l Vf the results for the phase shifts in the 3 D\ channel were close to 
zero, S( 3 Di) s» 0. The 3 Pq channel is displayed to exemplify that the P-wave interactions of the potentials are indeed 
small. 

Finally, observe in table|TT]that our cutoff values do not constrain the effective range to unphysical values, as alluded 
to in the discussion of the Wigner bound in sect. IPTl 

Tableland fig.Q]demonstrate: (i) The quality of the fit is sufficient for a NLO calculation with the expected accuracy 
level of about 10% for EFTjf at NLO. (ii) The potentials in the 2-nucleon sector are approximately equivalent. 



V. RESULTS FOR THREE AND FOUR NUCLEONS 



We now present the results of our feasibility study of using the RRGM with EFTs, exemplified at NLO calculations 
in EFT 4. The section is divided into four subsections. The first reports on an analysis of the correlation between 
the triton charge radius and its binding energy. The second looks at the splitting of the trinucleon binding energies 
due to Coulomb effects. The third discusses the Tjon correlation, and the fourth contains the findings for the 
3-Helium-neutron scattering system. 



A. Three nucleons: triton charge radius 

Friar et al. [46[ observed a correlation between the triton binding energy B{t) and the triton charge radius using 
various nuclear force models, before Platter et al. [41] showed by a variation of the three-body interaction parameter 
at leading-order that it can be understood as a consequence of universality in EFT^. With the model spaces defined 
in sect. IIII1 B(t) and the corresponding ground state wave function (f\t) were calculated at LO and NLO, and from 
that the root mean square charge radius r* h of the triton: 

^=((i|E^1(l + rf)|^ 2 , (13) 

where r*j is the position and (1 + t 3 ) the charge operator of the i-th nucleon. The results in fig. [2] confirm the expected 
behavior of an increasing r* h for more loosely bound systems. Results of two EFT^ LO calculations are shown. The 
RRGM LO potentials map out a band (gray shaded area) which includes the datum. The lower and upper LO line 
from [41] (dashed gray lines in fig. result from a fit of the LO LECs C\® to either a Sl a s or B(d), a s , respectively, 
providing some measure of higher order effects. For the RRGM LO calculation, we used a St t to determine the LECs, 
and varied the cutoff from 200 MeV (top edge) to 1.6 GeV (lower edge) in steps of 100 MeV. We choose the lower 
bound for A to be in the region of the pion mass. The upper bound is set by the observation that there is essentially 
no change for cutoffs larger than 1.6 GeV. The band appears to be saturated, with the bulk of its width coming from 
the region A € [400; 800] MeV. For fixed A, the 3NI was finally tuned to a given B(t) to generate the correlation lines. 

The quasi-exact Faddeev calculation of [Jl[ and the RRGM results do not overlap. Since app. O demon- 
strates that numerical inaccuracies of the RRGM are negligible, we speculate that this can be traced to the dif- 
ferences between the regularization schemes employed. The Faddeev calculation uses the separable cutoff function 
fniPyP 1 ) = exp (— -p 1 1 kF) exp (— p* 2 /k 2 ) with k > 1600 MeV, while the RRGM uses the non-separable regulator 

fh.{Pii?) — exp (p — p 1 ) 2 / A 2 ^j , plus an implicit regulator imposed by the finite number of width parameters. 

Therefore, both methods combined can be viewed as providing a check of residual regularization-scheme dependence. 
A conservative estimate of LO effects is thus the range of results coated by the combination of both methods. The LO 
accuracy at the physical triton binding energy is thus ±0.6 fm, and the measured charge radius happens to lie right 
in the middle of the LO band. Additional regularization schemes, e.g., with a range of cutoff values in the Faddeev 
approach, are thus conjectured to lead to an overlap between the results of [4l| and the RRGM bands. 
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B(t) [MeV] 

FIG. 2: The correlation between the triton charge radius and its binding energy. The narrower band mapped out by the 
EFTjf NLO results from the RRGM (blue squares, black dashed lines) as described in the text, compared to our RRGM LO 
results (gray shaded area) and a LO calculation by [4l[ (dashed gray lines), demonstrates convergence, also to the datum (r* h 
from [42§, B exp (t) from Q). The values from AV18(+ULX) (black triangles, RRGM calculation for this work) and a variety 
of other potential models (gray diamonds [3, 53) are consistent with our proposed EFTjf NLO correlation band. 

At NLO, the potentials without 3NI map out a correlation band which is more narrow than its LO counterpart, 
and again contains the datum. The width of the band can be estimated by considering the difference between the 
results using the extreme cases provided by potentials 2 Vf and & Vf. They produce the same NN scattering lengths 
and nearly the same triton binding energy, and their deuteron binding energies differ by less than 10%, while their 
triton charge radii differ by about 10%, consistent with the expectation of a NLO calculation. Finally, they are based 
on two significantly different cutoff values, A( 2 V^) « 400 MeV and A( 6 V^) ss 650 MeV. We therefore can base an 
estimate of the NLO band on the range mapped out by varying the 3NI for these two potentials. These lines are 
included in fig. [2] 

Variation of the cutoff, of the fitting-input, and of the 3NI lead therefore all to similar assessments of the uncertainty 
of the theory at NLO. At fixed B(t), the charge radius varies from LO to NLO by < 30%, in agreement with the 
power-counting which estimates the correction from order to order to scale as Q ~ ptyp/Aj, ~ ^. 

All three values, namely, ±0.2 fm from the nai've estimate Q « ^ and the observed convergence from LO to NLO, 
and ±0.1 fm from the above mentioned difference between 2 Vf and 6 Vjt, would be equally valid estimates for the 
theoretical uncertainty of this NLO calculation. As demonstrated in app.^l errors induced by numerical inaccuracies 
are negligible. Using the more conservative error estimate and the assumption for the NLO error band center at 
1.6 fm at the experimental B(t), EFT^ predicts a value 

r*f LO = (1.6 ± 0.2) fm (14) 

within the band of the leading-order value r 1 ^ = (2.1 ± 0.6) fm as quoted from [4l[. The NLO value is found in 
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good agreement with experiment [Hj], r^ p = (1.63 ± 0.03) fm. 

Another argument in favor of our definition of the correlation band is provided by the results of the two phe- 
nomenological models, AV18 and AV18+UIX [47]. In general, EFT predicts that the results of a potential which 
reproduces or shares input observables at least to the accuracy required at the considered order, deviate from the 
results of an appropriate EFT potential by less than the theoretical uncertainty of the EFT values at this order in 
the applicability range of the EFT. This criterion is easily met by AV18(+UIX), and hence its predictions have to be 
consistent with the proposed correlation band. The two-body potential AV18, reproducing the Nijmegen phases much 
more accurately than required to fall into this category of potentials, is expected to yield a value at a position within 
the {r l ch — B(t))-b&nd. The prediction for the triton charge radius of AV18+UIX, however, is expected to deviate less 
than 10% from the experimental datum within the error band, because this model has a three-body interaction added 
to reproduce the experimental B(t). Both expectations are consistent with the results shown in fig. [3J 

In conclusion, the results for these 3-nucleon observables show that although the potentials are approximately NN 
phase shift equivalent, they differ in their predictions of three-body observables. As mentioned in the introduction, a 
proper renormalization of the theory requires therefore one 3-nucleon contact interaction. Setting this 3NI to zero in 
the potentials 1-9 V^f , the expected dependency of observables in A > 2 systems on how the unobservable short-distance 
physics is modeled is observed. Different short-distance physics is modeled by the potentials not only by varying cutoff 
values but also by differing sets of LECs for the same cutoff, while the scheme- and regulator-dependent three-body 
interaction parameter is chosen as zero. With this parameter fitted to the triton binding energy, the prediction for the 
triton charge radius is consistent with experiment within the expected uncertainty range. A significant convergence 
from LO to NLO is observed. 



B. Three nucleons: effect of the Coulomb force 

In fig. [3l the RRGM results for the splitting between the binding energies of 3-Helium and the triton are shown 
for the potentials 1_9 V^. At NLO in the pion-less EFT, the strong interaction is isoscalar and hence does not break 
charge symmetry. Therefore, charge symmetry breaking (CSB) comes in our RRGM calculation only from including 
the Coulomb interaction between the protons in 3 He. The EFTjf results, hence, only show the model independent 
contribution of Coulomb interactions to the trinucleon binding energy splitting. The correlation band is mapped out 
by the EFTy. potentials with zero 3NI and by smoothly varying this three-body parameter for potentials > Vf. Both 
approaches result in correlations consistent with each other and with the results of the other NLO potentials. At the 
experimental triton binding energy, this leads to predicting 

(B(t) - S( 3 He)) = (0.66 ± 0.03) MeV . (15) 

Here, we estimate the theoretical uncertainty by the spread of the phase-equivalent NN potentials as in the previous 
section, e.g. comparing 3,9 Vf. The a priori error estimate at NLO of < 10% gives a larger uncertainty of ±0.07 MeV. 
However, one should keep in mind that including iso-scalar strong interactions by higher order terms of the effective- 
range expansion has identical effects on the strong interactions inside the triton and 3 He. These effects cancel out in 
the difference and only survive indirectly, as the strength of the Coulomb interaction in a system is also correlated to 
its size. We therefore quote the width of the correlation band as error estimate of our calculation. 

This value deviates by about 0.1 MeV from the experimental value of 0.764 MeV, see e.g. the recent review on CSB 
and Charge Independence Breaking (CIB) [50j. In line with the argumentation above, we attribute this difference to 
isospin breaking CIB/CSB interactions coming from the explicitly broken chiral symmetry in the strong and electro- 
weak sector from which only the parts resulting in the Coulomb force have been considered in this calculation. They 
enter in EFT^ only at higher order. To support this assertion, results of the potential models AV18(+UIX), which 
contain CSB interactions, are included in fig. [3] Both potentials are not elements of the correlation band suggested 
by the EFTjf points but agree with a shifted band, centered around the datum. In contrast, the values from charge 
symmetric potential models in fig. [3] lie within the NLO EFTjf band. This leads us to predict a model independent 
CSB/CIB contribution to the binding difference in NLO EFTjf at the experimental triton binding energy of 

{B{t) - S( 3 He)) CSB/CIB = (0.10 T 0.03) MeV , (16) 

anti-correlated with the Coulomb contribution to give the experimentally established difference. This is to be compared 
with the contributions from 2- and 3-nucleon CSB interactions which stem from Chiral Effective Field Theory, Breit 
and vacuum polarization corrections, and from corrections to the kinetic energy operator. In Ref. (49[, these were 
calculated to sum up to (0.112=f0.022) MeV, leaving about (0.652±0.022) MeV for the soft photon effects, dominated 
by the Coulomb interaction. This is in perfect agreement with the EFTy- result. 
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FIG. 3: Binding energy difference between 3 He and 3 H for EFT,f potentials (filled squares) compared to various isospin invariant 
potential models calculations [45l. |48| (circles) and RRGM values for AV18(+UIX) (triangles) which contain charge symmetry 
breaking terms. For the potentials 1 Vf (solid) and 7 V 1 t (dashed), a smooth variation of the 3NI leads to the two correlation 
lines. The upright cross is the experimental value without the contribution from CSB terms as found in [49l ]. 



C. Four nucleons: bound state 

Venturing into the four-body system, we now present results of the Tjon correlation line[l2| including Coulomb 
interactions between the ground state energies of the triton and 4-Helium at LO and NLO, extending the LO analysis 
officii], where the effect of Coulomb interactions was only estimated. 

Before reporting the results, we describe the employed variational space and its construction. The model space in 
which the 4-Helium binding energy B(a) was initially calculated is spanned by vectors defined in eq. (fSJ with relative 
angular momenta Iki < 2 including all coupling schemes to yield a total angular momentum state J 7 ' — + and width 
parameters "/dk to allow for the formation of triton, 3-Helium, and deuteron fragments. As mentioned above, the wq3 
sets were used for the deuteron, resulting in 169 configurations. For the 3-nucleon fragments, a small model space 
of dimension d < 70 was generated for each potential as follows. To an initial set of 20 basis vectors which bind the 
triton, one new basis vector was added, and its two width parameters ^dk (see eq. ([3])) optimized with the genetic 
algorithm to maximize the gain in binding energy. This process was iterated until \B 224 (t) — B small (t)\ < 500 keV, 
where the superscript labels the aforementioned 224-dimensional model space. For 3-Helium, a copy of this model 
space was used which differed only in its isospin quantum numbers. On average, 300 configurations were considered 
for the model space of a given Vf , for each of which four or five inter- fragment width parameters were used, resulting 
in roughly 1200 basis vectors. 

To assess if the reduction in the number of triton components still guaranteed for an almost complete 4-Helium 
model space, the binding energy B(a) thus obtained was compared to one calculated in the much larger scattering 
model space. Nevertheless, the two values differed by less than 50 keV. This model space, used for the + — channel 
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FIG. 4: Correlation between the triton and 4 He binding energies (Tjon line). The filled squares are the results using NLO 
EFT^ NN potentials from table U with no 3NI, and the empty squares represent the predictions for the 3NI fitted to B(t). The 
shaded NLO band results from a variation of the 3NI for those potentials. We compare to LO results from [lOfl. where the 
upper (lower) dashed line was calculated with LEC fitted to a s ,t (B(d),a s ), and to our RRGM calculation for A ranging from 
400 MeV (lower edge) to 1.6 GeV (top edge of light gray area). B cxp (a) is taken from [5l[, the values for AV18(+UIX) are 
reported in Jg^], and those using a chiral expansion to NLO and N 2 LO in 0. 



scattering calculation is denned below and is spanned by more than 7000 basis vectors. This study of reducing the 
dimension of the 3-nucleon fragments in the a-particle without deviating significantly from the assumed converged 
value is crucial for future applications of the method to A > 4 systems. 

In fig. 21 the results are compared to the LO band calculated by Platter et al. [10]. The spread of the NLO values 
is not in conflict with EFTj which allows a < 10% uncertainty at NLO. As for the triton charge radius, Platter et 
al. obtained the upper (lower) boundary of the LO band by choosing different NN observables, a S7 t(B(d) 7 a s ), to fit 
the LO LECs for cutoffs high enough so that B(a) did not change when further changing the cutoff. They already 
pointed out that this gives only a crude estimate of higher order effects. Figures [2] and [4] confirm that the LO accuracy 
would be overestimated with this method. Our LO calculation uses the same interactions as in sect. IV Al with a S: t 
as input, and cutoff values from 400 MeV (lower bound of LO band in fig. 2]) to 1.6 GeV (upper bound). In contrast 
to our RRGM LO correlation between B(t) and the triton charge radius, the width of the RRGM LO Tjon band is 
however not converged. The shifts between the positions of the correlation lines corresponding to cutoffs of 400 MeV 
(lower edge of gray area in fig. @|, 700 MeV (dashed dotted line), and 1.6 GeV indicate that a variation of A even 
beyond 1.6 GeV would be necessary to assess the LO uncertainty from cutoff variations only. However, elaborate 
technical modifications are required for the RRGM calculation at higher cutoffs. Analogous to the conservative 
estimate for the LO uncertainty in fig. [5J the LO Tjon correlation band is thus mapped out by both the RRGM and 
the Faddeev/Yakubovsky (thick dashed lines) results, which overlap nicely in this plot. This combined correlation 
band includes the datum and the narrower NLO band, which, from fig. 21 has at the experimental B(t) a width of 
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about 5 MeV centered around 28 MeV and results in a prediction of 

B NLO (a) = (28 ± 2.5) MeV , (17) 

which is consistent with the expected NLO uncertainty of about 10% and with experiment. Again, the results of the 
AV18(+UIX) models lie within the proposed band as it is expected of all interaction models of at least NLO. The 
observed broadening of the correlation band is a manifestation of the momentum dependence of the EFT expansion. 
The accuracy decreases with increasing typical momentum, eventually leading to a breakdown of the expansion. 

From the fact that there is still a one-parameter correlation, we conclude that no four-nucleon contact interaction 
is required to renormalize the theory at NLO. One three-body parameter fitted to data suffices to yield proper NLO 
predictions for four-body observables within the theoretical accuracy (empty squares in fig.|4]). 



D. Four nucleons: scattering 



We now turn to scattering observables. In principle, all low-energy observables should be correlated with the triton 
binding energy. The recent results [53[ for the singlet and triplet n— 3 H scattering lengths using three potential models 
and a N 3 LO chiral potential are evidence for this assertion in the four-nucleon scattering system. Here, the real part of 
the S-wave spin singlet scattering length ao( 3 He-n) for elastic 3-Helium- neutron scattering is investigated. In fig.[5l its 
value is shown as a function of B{t) for six potentials V* including the Coulomb interaction. For two NLO potentials 
with A = 440 MeV (solid line) and A = 550 MeV(dashed line), with 2-nucleon LECs fixed, we also show the effect of 
a smooth variation of the 3-nucleon interaction parameter. 

To extract ao( 3 He-n), six two-fragment channels, 3 He-n, t-p, d-d (Z re i = 0, 2), (nn)-(pp), and dq-dq (dq is the singlet 
"deuteron" with 5 = 0), are included. The latter two consist of unbound fragments and usually model possible three- 
and four-body breakup reactions. They, as well as the two d-d channels, are however for this calculation only needed 
to provide configurations for distortion channels since only the t-p channel is open a few eV above the 3 He-n threshold. 
For the fragment wave functions tpj (see eq. (|4])), a 224 dimensional basis was used for the triton and 3-Helium, and a 9 
dimensional one for the deuteron, whose six Lj = vectors built the nn, pp, and dq states. For these six channels, the 
20 width parameters W12 were used for the ojj m in eq. Almost all configurations included to build those physical 
channels could be recycled as distortion channels to allow for more freedom in the minimization of the variational 
functional. Less than ten configurations had to be excluded to avoid numerical linear dependences. In each distortion 
channel, four to six relative width parameters u>j m , taken from W12 with u)j m > 0.02 fm -2 , were used. Numerical 
stability and convergence of cto( 3 He-n) were assessed by increasing the number of included relative widths Uj m by one 
for each distortion channel, yielding changes in ao( 3 He-n) of the order of the numerical uncertainties, given that the 
initial u)j m were chosen appropriately. The lowest eigenvalue of the Hamiltonian is equal to the ground state energy 
of 4-Helium in this model space and was allowed to change in this process by not more than 10 keV. Significantly 
larger changes in this eigenvalue which leads to a result not of the order of magnitude suggested by the LO Tjon band 
signal numerical linear dependences. If the model space is too small, or if the width parameters for the relative wave 
function were chosen inappropriately, changes of the order of 100 keV up to a few MeV are expected. 

The 3- Helium- neutron scattering length ao( 3 He-n) was calculated from the corresponding diagonal element S22 of 
the S-matrix at a center of mass energy of E cm < 10 eV above the 3 He-n threshold analogously to [541 ]: 

ao(E cm )= , 1 - S22 ^"°) . (18) 

/ fA/£ cm (l + S 22 (£ cm )) 

A fraction of the flux is diverted into the open triton-proton channel, resulting in a nonzero S12 S-matrix element 
and hence a nonzero imaginary part of ao( 3 He-n) which is not reported here. The 3NI variation yields an almost 
linearly increasing imaginary part with increasing B(t) from about Im{ao( 3 He-n)} « —7.5 fm at B(t) « 5.5 MeV to 
Im{ao( 3 He-n)} ss —2.0 fm at B(t) sa 9.1 MeV. This qualitative observation is consistent with the decreasing trinucleon 
binding energy splitting with decreasing triton binding energy (see sect. IV Bp . which results in a smaller separation 
between the respective thresholds in four-nucleon scattering. This handle on the threshold separation suggests an 
approach to circumvent numerical problems associated with the proximity of thresholds by extrapolating results for 
their physical values from calculations performed at more deeply bound 3-nucleon states. 

In fig. [5j the predictions of the NLO potentials with and without 3NI for Re{ao( 3 He-n)} decrease with increasing 
triton binding energy. They map out a band which includes the datum. A comparison to hard sphere scattering 
qualitatively explains this behavior. A higher 3-Helium binding energy i?( 3 He) corresponds to a smaller nucleus 
analogous to the triton as shown in fig. [2j As the scattering length is proportional to the radius of the hard sphere, 
a ( 3 He-n) is expected to decrease for increasing _B( 3 He). We define the NLO correlation band to be centered around 
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FIG. 5: The correlation between the triton binding energy and the real part of the spin singlet, S-wave scattering length 
ao( 3 He-n) for elastic 3 He-n scattering. The dots represent the values of the NLO EFT^ potentials with a numerical uncertainty 
of 0(0.02 fm). A variation of the 3NI yields the continuous correlation lines for two selected NLO potentials with a fixed cutoff. 



the values predicted by the potentials with zero 3NI. The change in ao( 3 He-n) observed between 2 Vf and 6 Vf serves 
as an error estimate for the band centered around 7.5 fm. Again, the assumption of the center of the band at the 
experimental B(t) follows from the values in fig. [5] Thus, EFT^ LO reports 



Rc{a ( 3 He-n)} = (7.5 ± 0.6) fm 



(19) 



The 3NI was varied only over the depicted range of trition binding energies, where one and the same fixed RRGM 
variational space can be used. The fragment model space for the two interactions was optimized for the 3NI parameter 
fitted to B(t). Only the number of included distortion channels was changed to reach convergence within this space for 
each value of the 3NI. The band mapped out by the two 3NI lines for (7.5 MeV < B(t) < 9 MeV) includes the datum, 
and its slope is consistent with the one indicated by the 6 NLO potentials without 3NI. This explicitly demonstrates 
that a variation of the 3NI has the same effect as varying the short-distance part of the NN interaction. The NN 
model AV18 yields a value within the error band. Adding the UIX three-body interaction moves this point into the 
10% NLO uncertainty radius around the datum. The numerical accuracy of this scattering calculation is assessed in 
app.|B]to be better than 1%. 

The conclusion drawn from the newly found correlation in fig. [5] is that every potential with the correct NN low- 
energy phase shifts and appropriately tuned three-body interaction, e.g., to give the correct triton binding energy, 
predicts not only the correct B(a) but also the experimental ao( 3 He-n) within a NLO error range. In table [TTTl the 
results for a ( 3 He-n) of RRGM calculations with the EFT^ LO potential and the AV18 NN force model with and 
without the three-body interaction models Urbana9 (UIX) and Illinois2 (IL2) (55[ are given. 

The values for the recent IL2 3NI have been calculated for this work, employing methods described in [H| • The 
IL2 prediction for Re{ao( 3 He-n)} is almost identical to the UIX value and not plotted separately in fig. [5J UIX and 



17 



TABLE III: RRGM predictions for the triton binding energy and the spin singlet, S-wave scattering length ao( 3 He-n) for elastic 
3 He-n scattering of phenomenological nucleon potential models. The values for AV18 and AV18+UIX are taken from [54 |. while 
the AV18+IL2 numbers are new results of this work. The imaginary part was investigated on a qualitative level, only. 



force 


B{t) [MeV] 


Re{a } [fm] 


Im{a } [fm] 


eft^ u 


8.48 (input) 


7.5(6) 


-2.6(?) 


AV18 


7.57 


7.81 


-4.96 


AV18+UIX 


8.43 


7.62 


-4.07 


AV18+IL2 


8.48 


7.63 


-4.28 


exp 


8.48 


7.456(20) |56J 
7.370(58)^3 





IL2 have parameters fitted to, amongst others, the triton and 4-Helium binding energies. The observed deviation of 
the prediction of both models for ao( 3 He-n) is therefore not a result of a deficiency of the structure of the potential. 
The correlation between B(t) and ao( 3 He-n) supports the conjecture that this small deviation can only be improved 
by the inclusion of higher order interactions. 

Analogously to UIX, the IL2 model easily satisfies the criteria mentioned in sect. IV A[ namely to reproduce low- 
energy NN observables and the triton binding energy at least with a 10% accuracy, and therefore should with its 
prediction for ao( 3 He-n) also lie within a 10% radius of the datum in fig. [5] This EFT prediction is confirmed here by 
explicit calculation. 

Both experimental values for eto( 3 He-n) are included in the predicted universality band, and therefore we cannot 
resolve the discrepancy between the two measurements. 

E. Evidence for the absence of four-nucleon interactions at NLO 

We close with a comment on the implications of these results for the scaling of four-nucleon interactions (4NI) in 
EFTy- at LO and NLO. 

Simplistic dimensional analysis suggests that a momentum-independent 4NI enters at N 4 LO, but the unusual 
renormalization of the 3NI may also promote the 4NI to contribute at lower orders to ensure renormalizability. By 
varying their harsh cutoff between 800 MeV and 2000 MeV, Platter et al. [10( concluded that a 4NI is not necessary 
at LO for cutoff-independence of the 4-Helium binding energy. 

In the present approach with a Gaussian regulator, the cutoff was varied for the 4-Helium binding energy and for 
the 3-Helium-n scattering length from 400 MeV to 700 MeV. This led to a correlation band near-identical with the 
ones mapped out when using different NN potentials (in part also at different cutoffs) or different 3NIs. All results 
were converged numerically. Moreover, we demonstrated that the theoretical uncertainty for B(a) decreases from LO 
to NLO by a factor consistent with the a-priori expansion parameter estimate Q ss i. In both cases, the physical 
datum lies well inside the NLO correlation band. 

Let us assume now that a 4NI enters at NLO. Our not accounting for it would then have two effects: The results 
could be unstable against cutoff variations, indicating that a 4NI is necessary at NLO to renormalize EFTj. We see 
no such effect, but are aware that the cutoff might have to be varied beyond the window chosen in this exploratory 
study. Secondly, the theoretical accuracy of the results would be reduced to LO, even if results independent of the 
cutoff can be achieved, i.e., even if the cutoff can be removed to infinity without a 4NI. In other words, when one 
separately or combinedly varies the cutoff or the 2-nucleon potential used or the 3NI strength, the different results 
should spread in a corridor set by the size of LO corrections, i.e., < 30%, and not by the corridor of < 10% expected 
in a NLO calculation. We have demonstrated above for both observables that the residual short-distance dependence 
of the observables is < 10% of the central value, taking as conservative estimate the combination of errors which 
occur when varying EFT^ at unphysically short distances: Different, phase equivalent NN potentials; different 3NIs; 
different cutoffs. The correlation bands and error estimates are thus on the quantitative level consistent with those 
of a NLO calculation. 

We also note in passing that the physical datum lies well inside the NLO correlation band for both four-nucleon 
observables. A residual 4NI not necessary for renormalization but with unnaturally large coefficient can therefore be 
ruled out. 

Our conclusion therefore is: There is strong evidence that 4NIs do not enter at either LO or NLO in EFT^. An 
indisputable criterion to assess whether our interpretation is correct can be provided in a future study [5Sjj of the 
convergence pattern of the momentum-dependence of four-nucleon observables from zero to the breakdown scale. In 
it, an observable calculated up to order Q n must show a residual short-distance dependence on compatible with n 
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powers of the typical low-energy momentum. This analysis will not re-sum the effective range contributions into the 
2-nuclcon propagator, but treat higher order effects in strict perturbation. 

VI. CONCLUSIONS 

The effective held theory formalism can be used to explain empirically found correlations amongst few-nucleon 
observables like the Phillips- and Tjon line. In these two cases, it relates the deviations from data to, first, an 
incomplete renormalization in the 3-nucleon sector, and second, to higher order interactions omitted in the course of 
the EFT expansion. The theoretical uncertainty at every order of the calculation can be quantified in this approach. 
Numerical inaccuracies were demonstrated to be negligible in app. Eland app.|B] Different results for different short- 
distance parameterizations therefore have other origins. The size of the expansion parameter pty P /Ai, determines how 
fast an EFT expansion converges and decides its usefulness for the calculation of an observable in a given system. 
In heavier systems like 4- Helium, the expansion parameter can a priori be as large as 1, so they are border line. 
However, we find that next-to-leading order corrections to leading-order results are still parametrically small in the 
four-nucleon system. This confirms a pattern already seen in 2- and 3-nucleo n sy stems at momenta which approach 
the a priori breakdown scale but where convergence is still found, see e.g. 0, E3, Hl| ■ 

The computational challenges of calculating few-nucleon observables are met in this work by the Refined Resonating 
Group Method. It provides a versatile method for bound- and scattering properties. Here, it was found highly 
economical with respect to computer time when combined with the EFTy NN potentials derived in sect. IIVI for a 
range of momentum cutoff values and with a full treatment of the Coulomb interaction. 

For two correlations for which leading-order calculations exist, namely between the binding energy and charge 
radius of the triton in sect. IV Al and the binding energies of 4-Helium and the triton in sect. IV Cl our coordinate space 
EFTy. calculations at NLO report the expected improvement from LO to NLO consistent with an expansion parameter 
Ptyp/Ab s=s |. By that we demonstrated that a consistent description of the a- particle is possible at NLO in EFTy . In 
sect. IV 51 we also report a correlation between the triton binding energy and its difference to the 3 He binding energy. 
As the EFTjf potential is at NLO iso-spin symmetric, this model independent difference is attributed to Coulomb 
interactions only, which are included in the RRGM. At the physical triton binding energy, this value agrees well both 
in magnitude and uncertainty with estimates of charge-symmetry breaking and Coulomb contributions to 3 He binding. 
In sect. IV Dl a new correlation between the triton binding energy and the real part of the singlet S-wave scattering 
length of 3-Helium-neutron scattering similar to the Tjon line is also found. This, and the three aforementioned 
correlation bands, let us also conclude that no four-body contact interaction is required to renormalize the system at 
next-to-leading order. The position of all four bands, which represent universal properties of the 2-nucleon system, 
was determined by fitting nine NN potentials differing at short-distances but with identical long-distance behavior, 
by variations of the 3NI strengths, and by changing the cutoff. 

Consistent with a basic tenet of EFT, namely model independence, we also showed that the results of the phe- 
nomenological models AV18(+UIX/IL2), which share the input of our EFTy. NLO potentials, agree with their results 
within NLO accuracy. 

Future work will utilize the relatively fast computations of four-nucleon observables, which result from the harmonic 
interplay of the RRGM with EFT^ LO potentials, in analyses of universal properties of A > 4 systems. The Borromean 
halo nucleus, 6-Helium,_is of special interest because EFTy. calculations in this system can also provide input for aN 
effective field theories [59j . Furthermore, we are now equipped for a study of electro- weak interactions with heavier 
nuclei like 4 He(7,p) and 4 He(7,n) using EFTy. and the RRGM, to address both conflicting measurements as well 
as theoretical calculations [60- 62] . Valuable input for astrophysical calculations for the prediction of light element 
abundances can also be provided by a calculation of reaction cross-sections for e.g. d(d,n) 3 He and d(d,p)t at very low 
energies. 
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Appendix A: Numerical stability: two-nucleon sector 

In the analysis of sect. [V] various methods are used to ascertain higher order effects and provide reliable theoretical 
uncertainties of our NLO calculations. It is therefore imperative to ensure that purely numerical inaccuracies of the 
variational method are not a significant source of error. 

When increasing the model space, two obstacles have to be considered: first, numerical linear dependences amongst 
the basis vectors; and second, too broad distortion channels. Both issues lead to unstable results and occur especially 
when the dimension of the model space becomes large. The bulk behavior of the phase shifts in the 4-Helium system 
was stable in that respect, but the prediction for the 3 He-n spin-zero S-wave scattering length was much more sensitive 
to the addition of basis vectors. This is a consequence of too broad distortion channels, which is associated with the 
expansion of the Coulomb functions. We illustrate this point here as a purely numerical issue, anticipating results 
from sect. |Vj All integrals needed to calculate a\j and b\j m have a support only for Rj < i? ma x, where i? ma x is set 
by the size of the largest fragment. When the Coulomb functions are expanded in Gaussians, this inner region is 
weighted more at the expense of a less accurate fit for larger separations. This poses no problem if broad Gaussians 
are used for the distortion channels, i.e., uij m < w m i„(i? max ), because the results are only affected by the values of 
the Coulomb function in the outer region, Rj ^> i? max . In sect. |V1 potentials are presented which yield relatively 
weakly bound 3 H nuclei with more extended wave functions (see fig. [2j. For those interactions, distortion channels 
with broader width parameters had to be included in order to form the triton cluster inside of the 4-Helium nucleus 
and reach convergence for the 4-Helium binding energy. If, on the other hand, too many of those broad vectors are 
included, the phase shifts and especially the scattering length ao( 3 He-n) become unstable. 

The NN model space used to fit the LECs (see sect. IIV[) is a key factor. If the scattering and bound states of 
a given set of parameters cannot be expanded accurately, the resultant potential will carry a non-negligible model 
space dependence. To minimize this dependence, a large two-body model space was chosen with a set of width 
parameters which cover an interval corresponding to the expected extention of the 2-nucleon wave function. Twenty 
width parameters in the range ujj m € [0.001 fm -2 ; 130 fm -2 ] (called set W120 [sH) were used for the NN scattering 
states, and the same set for the widths jdk in eq. © for the S- and D-wave components of the deuteron bound state. 
Completeness of this set was defined to be sufficient by comparing the results (see table [TTJ) to those of another set 
of 20 parameters, where each inverse width of W120 is divided by ten, resulting in broader widths. Hence, all widths 
of this set W12 lie in the interval [0.0001 fm -2 ; 13 fm -2 ]. For the triton and 3-Helium wave functions, a previously 
mentioned set of 224 vectors taken from [54| spanned the model space. It contains all possible couplings with total 
angular momentum L < 2 and 142 different width parameters. To model the deuteron fragment in the 4-Helium 
calculation, its model space was reduced to 6 S-wave and 3 D-wave widths for each potential separately, employing a 
genetic algorithm which optimizes the width parameters of the basis vectors [34( . These sets are named wq3 in what 
follows. These optimized deuteron width sets for narrower potentials with larger cutoff values also contained narrower 
widths, in line with the need to model wave functions with a more complicated short-distance structure. 

The expansion of the wave function in terms of Gaussians had significant impact on the choice of the regulator. 
First, a cutoff momentum space regulator f(q) = exp (— q 1 /A 2 ) , leads to a Gaussian radial dependence as shown below 
in eq. (|5]). Hence, all radial dependences of the potential are Gaussians and can directly be implemented into the 
RRGM. Second, with its functional form fixed, the magnitude of A fully specifies the regulator and is also influenced 
by the RRGM. A lower bound is set by the breakdown scale A& « of EFTj, since modes which lie in the range of 
applicability of EFT^ should not be suppressed. On the other hand, an upper bound A t for the cutoff is set by the 
numerics of the RRGM. Large cutoff values correspond to narrower Gaussian potentials in coordinate space, and hence 
the RRGM model space must include narrower Gaussians as well to expand the corresponding, more localized NN 
bound state accurately, see below. Clearly, A cannot be chosen considerably larger than the inverse of the narrowest 
width in which one expands. A LO calculation in the deuteron channel, i.e. V(r) — Ctcxp (— A 2 f 2 /4), was carried 
out to illustrate this point. The low-energy constant Ct was fitted to reproduce the experimental neutron-proton 
triplet scattering length <Z( Xp by numerically solving the two-body Schrodinger equation with MATHEMATICA. The 
scattering length was calculated from the phase shift at E cm — 10~ 4 MeV, and the resulting coefficient Ct was then 
fed into the RRGM code. The RRGM result for the deuteron binding energy is plotted against A in fig. [5] The 
LO prediction of the effective range expansion, B(d) = 1.4096 MeV from o,( Xp = 0.02748 MeV -1 38], is reproduced 
in both considered model spaces consistently in the RRGM calculation up to the cutoff A t . Above this threshold, 
AJ" 12 ps 1.4 GeV and AJ" 120 ps 3.6 GeV, the RRGM solution becomes strongly cutoff-dependent because the respective 
model space is insufficient to expand a bound state corresponding to such narrow potentials. We explain this as 



20 



2.5 



„ 1.5 

> 
CD 



CO 1 



0.5 







model space 1 + 
model space 2 X 
B d (ERE) = (a t 2 M)" 1 — 




1000 



2000 3000 

A [MeV] 



4000 



5000 



FIG. 6: Binding energy Bd of the deuteron in leading-order EFT^ as a function of the Gaussian momentum cutoff A for two 
different RRGM model spaces. The space 10120 contains narrower width parameters, and hence can be used for the potentials 
corresponding to larger values of A. Using the triplet neutron-proton scattering length as experimental input, the effective 
range formula predicts Bd — 1.41 MeV. 



follows. The radial ground state wave function reaches a maximum around r p w A -1 , before decaying exponentially 
outside the interaction region. Most economically, this peak can be approximated by a sum of two Gaussians with 
widths comparable to r p . For the expansion of the exponential decay, broader Gaussians suffice. If the model space 
does not contain Gaussians of width r p , an approximation of the increasingly steep rise of the wave function at zero 
accompanying an increasing A will eventually fail. Hence, bound states can no longer be expanded in the model 
space for potentials of shorter range than the narrowest Gaussian basis state. In the case shown in fig. [51 relating 
the narrowest Gaussian widths w max of the two model spaces, and Wi2Q, via the regulator exp (— A 2 r^/4) to a 
cutoff A t = 2^/w m ax, one expects the ground state in the respective model space to become unbound for cutoffs larger 
than A™ 12 ss 1.4 GeV and A™ 120 4.5 GeV. The thresholds suggested in fig. [5] by a drop of the deuteron binding 
energy are of this expected magnitude. This lead to the conjecture that the W120 RRGM model space is appropriate 
for calculations with cutoff values in a range 150 MeV < A < 3 GeV. 

A comparison of the phase shifts 5( 3 Si) below E cm — 10 MeV resulting from an RRGM calculation on the one side, 
and of a numerical integration of the Schrodinger equation on the other strongly supports this. The results shown 
in fig. [7] of a potential with a cutoff of 1.5 GeV differ by less than 0.3%. As long as the cutoff was kept below the 
threshold At, the relative difference was always found to be of that order. For cutoff values approaching or surpassing 
At, the discrepancy increased considerably indicating the predicted failure of the specific RRGM model space. 

For the NLO potentials, A was taken from the interval between 400 MeV and 1 GeV, in which w\2 and W120 yielded 
the same deuteron binding energy. Therefore, it is reasonable to assume that both model spaces are large enough to 
expand the scattering and bound states of the NLO potentials, so that the model space dependence of the potentials 
should be minimal. Differences between results using different potentials and short-distance physics do therefore not 
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FIG. 7: Relative deviation between the Si NN phase shift results obtained with the RRGM and a numerical integration of 
the Schrodinger equation. A leading-order potential with one parameter fit to a t and a cutoff A = 1.5 GeV was used. Notice 
the scale on the ordinate. 



stem from numerical inaccuracies. 



Appendix B: Numerical stability: four-nucleon scattering 

To gauge the size of the error introduced by numerical uncertainties from diagonalizing the Hamiltonian for the 
calculations with zero 3NI, a ( 3 He -n) was calculated over a wider range of energies, E cra < 2 keV. Figure [5] shows 
Re{ao( 3 He-n)} as a function of the matching energy at which it is calculated for the potential 8 V^ . The fluctuations 
resulting in the band-like shape with a width of about 0.015 fm out of cto( 3 He-n) 9.415 fm are due to uncertainties 
associated with the diagonalization of numerically singular matrices. The error due to effective range corrections can 
be estimated by observing the difference, 

|Re{a (200 eV)} - Re{a (2 keV)}| « 0.01 fm . (Bl) 

We therefore conclude that the numerical precision of our result is about 1% and is considerably smaller than the 
theoretical accuracy of ~ 10% attributed to a NLO calculation. For AV18, the same analysis produced stable results 
with respect to numerical fluctuations, i.e. no visible band as in fig. [5] is found. This is attributed to the model space 
used for AV18, which was optimized for this potential and purged of states with large but mutually canceling overlap 
with the ground state. As mentioned at the beginning of this paragraph, this analysis applies to values obtained in 
an optimized model space only. 
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FIG. 8: The real part of the spin singlet, S-wave scattering length for elastic He-n scattering as a function of the center of 
mass energy above the 3 He-n threshold at which ao( 3 He-n) is calculated from eq. (|18|) . Notice the scale on the vertical axis. 
Taking the central value of the band as the actual value, its decrease from 200 eV to 2 keV is due to effective range corrections, 
while the band is a result of numerical noise. s Vf was used for this plot, with the other Vf potentials showing similar behavior 
with respect to the size of the higher order corrections and numerical fluctuations. 



Appendix C: Computing time requirements 

In tablc lIVl the computer time, normalized to one CPU, of the EFTjf and the AV18(+UIX) calculations is compared. 
We display the time used for the build-up of the Hamilton matrix, which surpasses by at least an order of magnitude 
the time needed for the other stages, namely the calculation of the spin- and coordinate space matrix elements and 
the diagonalization of the Hamilton matrix. The EFT potentials allow for faster calculations because we choose the 
regulator of the EFT potentials such that no expansion of the radial dependences is necessary. Hence, only one 
matrix element is evaluated for a specific operator. For the same operator which has a different radial dependence 
in AV18, one has to calculate as many matrix elements as Gaussian basis functions are necessary to approximate 
this dependence accurately. The calculational time per term in the Gaussian expansion is dominated by the operator 
structure only, and is found to be comparable for AV18 and the EFT potentials. This resembles the already stressed 
occurrence of all AV18 operators in the EFT potentials considered here. The central 3NI is equally inexpensive with 
respect to computation time for systems with A < 4 as the NN potentials. We plan to invest the computer time 
saved by the convenient choice of the regulator of the EFT^ interaction in calculating heavier nuclear systems [58j . 
For those systems, the model space and the number of possible spin and coordinate spaces coupling schemas increases 
the number of matrix elements to be summed up considerably. Having reduced the number of matrix elements from 
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TABLE IV: Runtime comparison between RRGM calculations with EFT and phenomenological potentials, dim(rad-dep) is 
the number of Gaussians needed to expand the radial dependences, and dim(ms) is proportional to the dimension of the model 
space. The parallel computations were performed on IA32 Xeon 2.66 GHz nodes of the RRZE in Erlangen, while for serial 
code, an E6750 architecture at 2.66 GHz was used. 



Observable 


Potential 


dim(r£icl-c 


ep) dirn(ms) 


Runtime 
[CPUh] 


Bit) 


AV18 


44 


67 
67 


0.04 


TTTYfrinrallpl^ 
kj iyv i pdi duel i 

EFT x V f 


9 


n 1 7 
< 0.01 




ll l/ / +3NI 


1 




< 0.01 


B(a) 


AVI 8 (parallel) 

UlX(parallcl) 

EFT 


44 
204 
9 


324 
550 


2.6 
6.8 
0.45 




1 V>+3NI 


1 




0.03 


a ( 3 He-n) 


AV18(parallel) 
UlX(parallel) 
EFT 1 V f 


44 
204 
9 


324 
550 


10.4 
27.1 
1.8 



a Gaussian expansion of the potential, therefore, allows for more refined and larger model spaces. 
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